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We study the quark number susceptibility of the hot quark-gluon plasma at zero and 
fT") ' non-zero quark number density, using lattice Monte Carlo simulations of an effective 

^vq theory of QCD, electrostatic QCD (EQCD). Analytic continuation is used to obtain 

results at non-zero quark chemical potential [i. We measure both flavor singlet (diagonal) 
and non-singlet (off-diagonal) quark number susceptibilities. The diagonal susceptibility 
approaches the perturbative result above ~ 20T C , but below that temperature we observe 
significant deviations. The results agree well with 4d lattice data down to temperatures 
^ ' ~ 2T C . The off-diagonal susceptibility is more prone to statistical and systematic errors, 

but the results are consistent with perturbation theory already at 10T C . 
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1. Introduction 



The quark (baryon) number susceptibility of hot QCD matter characterizes the "softness" 
of the equation of state. It is directly related to the event-by-event fluctuations observed 
in heavy ion collision experiments [T], probing the phase diagram and the properties of the 
hot QCD plasma. Thus, it is of significant interest to calculate it theoretically as accurately 
as possible. Hence, several calculations of susceptibility have been published using lattice 
simulations [HEHHEIEIEIIE] or perturbation theory [9iri0| fTT]. 

In this work we use lattice Monte Carlo simulations in order to measure the diagonal (flavor 
singlet) and off-diagonal (non-singlet) quark number susceptibilities at high temperatures and 
at non-zero densities. Instead of full 4-dimensional QCD, the theory we study on the lattice 
is a dimensionally reduced effective theory of the hot quark-gluon plasma phase of QCD, 
electrostatic QCD (EQCD) [HJ Q3J CGS QI] . It is by now well established that EQCD can 
accurately describe many properties of the hot QCD plasma, and it provides a very convenient 
starting point for studying high-temperature QCD using perturbative analysis |17l \18\ [19] or 
non-perturbative lattice simulations. 

The validity of the effective theory approach is based on the fact that at high enough 
temperatures the gauge coupling constant g becomes small, giving rise to three relevant 
momentum scales (neglecting quark masses): hard scale p ~ irT, corresponding to non-zero 
Matsubara frequencies, soft electric scale ~ gT and supersoft magnetic scale ~ g 2 T. EQCD 
is obtained by (formally) integrating over the hard scales perturbatively, leaving an effective 
theory for soft and supersoft scales. All infrared divergences inherent in finite temperature 
field theories are correctly contained in the effective theory. A crucial feature of EQCD is 
that all of the fermionic modes are integrated over, leaving a purely bosonic theory. 

EQCD offers an interesting alternative to standard high-temperature lattice simulations. 
Above all, the theory is three-dimensional and purely bosonic, making it much cheaper to 
simulate. The standard QCD lattice simulations work well at temperatures up to 5-10 
T c , but due to the sheer cost of the simulations with light quarks it can be very difficult to 
obtain accurate results. In contrast the perturbative analysis works at temperatures T > 10T C 
(albeit with slow convergence), but since the infrared singularities in the magnetic sector 
cannot be treated perturbatively the accuracy is limited to some order (depending on the 
observable) in the coupling constant expansion. The lattice simulations of EQCD fully include 
the effects of the infrared singularities, thus offering a clear way to improve on the perturbative 
results. While EQCD cannot describe the QCD phase transition, it has been observed to be 
quantitatively accurate down to temperatures of order 2-4 T c , depending on the quantity of 
interest. On the other hand, it is relatively easy to do EQCD simulations at arbitrarily high 
temperatures, enabling one to quantify the convergence to the perturbation theory and the 
role of the infrared singularities. Lattice simulations of EQCD have been used to calculate 
QCD pressure at high temperature [20\ I21j. spatial string tension |22j . and spatial screening 
lengths [231 121 EH]. 
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In this paper we present the lattice calculations using EQCD to measure the diagonal 
and off-diagonal quark number (baryon number) susceptibilities at zero and non-zero baryon 
chemical potential. At non-zero chemical potential EQCD suffers from a sign problem, al- 
beit this is milder than in full QCD. The finite chemical potential results are obtained by 
performing simulations with imaginary values of the chemical potential and then analytically 
continuing to real chemical potential. We observe that the deviations from the perturba- 
tion theory are significant up to temperatures of order 20T C . On the other hand, EQCD is 
observed to work at surprisingly low temperatures: our results agree well with existing 4d 
lattice simulations even slightly below 2T C . The method also is well suited for simulations at 
non-zero chemical potential, because our observations agree those of [5] and extend to even 
higher values of chemical potential. The results have been partly published in |26|, I27j. 

The paper is organised as follows. In Sec.[2]we give the theoretical background and specify 
the considered observables. In Sec. Owe present the numerical results of lattice Monte Carlo 
Simulations. Conclusions are given in Sec. [H 



2. Effective theory 
2.1. Action 

The electrostatic QCD with finite chemical potential fi is defined by the action 
Se = <1 3 x£e 



Ce = -Tr[4] + Tr[A,Ao] 2 + m|Tr[A2] + i 73 Tr[A|]] + A3(Tr[A2]) 2 , (2.1) 

where F^j = diAj — djAi + igs[Ai,Aj] and Di = di + igsAi. Fij, A{ and Aq are traceless 
3x3 Hermitean matrices (Aq = A^T a , etc). The theory has 4 parameters: g\ (3-dimensional 
gauge coupling), m|, A3 and 73, with dimensions [g 2 ] = [A3] = GeV, [73] = GeV 3//2 and 
1 = GeV 2 . Non -zero value of the parameter 73, caused by non-zero quark chemical 
potential, renders the action complex. Thus, this theory is not free from the sign problem of 
finite density QCD. 

It is convenient to define three dimensionless ratios 

y = —, x = —, z = -3, 2.2) 
93 9i 93 

leaving only g| dimensionful. Through the dimensional reduction process (perturbative 
matching of suitable observables in EQCD and real QCD), the parameters of EQCD be- 
come functions of physical 4d parameters: the temperature T and the chemical potential fi 
(the quark masses are set to zero) . The parameters are also functions of the renormalization 
scale Ajjjg used in the derivation of the effective theory. If we denote the number of quark 
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flavors by Nf, for N c = 3 the relations are |16|. I24|: 
9 24vr 2 T 

23 



i-(l-g^(W, + ^) (2.4) 

\A\n-.r \ 1 + ^6 + Wt^) + 

\ 1=1 / 

m-sm-uNj^/ * gjjfeg!) aft 'Vo M (2.5) 

96vr 2 (9-A f ) I ^486-33iV f -lliV f 2 -2iV f 3 ^M V 7 V ; 

^ fiw / 21 +3Af \ n N . , 



33-2iV f A s /A^ 

9 - iV f 
33 - 2iVf A^/A- 

(9-AT f )(6 + iV f ) X Ai 



3vr V 18 - 2A f 
where /2 = /j,/(ttT), and, for small /2, X>(/2) « — 7C(3)/2 2 /2, and 



i=l 

-,2 



- /-3 + 4A f log4 \ , . 

P V 594 — 75iVf + A 2 'J y ' 

The dimensional reduction scheme is expected to be valid temperatures down to ~2-4T c and 
chemical potential up to [i ~ irT or p, ~ 1. For JVf = 2 these values correspond to a; ~ 0.1 
and z ~ 0.1. At higher temperatures x becomes rapidly smaller. Hence, the higher order 
corrections in x in above formulas become in practice very small, and we ignore corrections 
0{x) in above expressions. We shall further restrict ourselves to 2 massless quarks, Af = 2: 

d = ffiU=o 

X = 2^=0 

y = V U=o ^ + I E Mil = 2/0 ^1 + I £ (2.9) 
^ 3vr 

See [24] for more discussion about the effect of this approximation. 
2.2. Susceptibility 

We define the quark number susceptibility in EQCD as: 

i d 2 , _ i a 2 



X3,ij 



V d/iidfij 



\nZ = ±--^—\n [ VA k VA exp(-S E ), (2.10) 
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where i, j stands for quark flavors u and d, and label 3 indicates that this is a result from 
3-dimensional effective theory. Thus, there are two independent components of the suscep- 
tibility: diagonal (i = j) and off-diagonal (i ^ j). Using the shorthand notation for the 
dimensionless volume averages 

A o = 4y f d 3 xTrA%(x), (2.11) 

and defining the condensates 

(h = (M) 

C2 = Vgt (<(ig) 2 > - (Al) 2 ) 

Cs = Vgi [{(Alf) - {Alf) (2.12) 
C, = V9t {{AlAl)-{Al){Al)) , 
we can write the susceptibility as 

^jiT = Vo Ci - g^2° 2 + J^Vifij Vo c 3 + ^(fii + N) Vo c ± (2.13) 

We note here the rather striking fact that the expectation value in C4 is purely imaginary 
for real p,, rendering the full expression real. The imaginary expectation value comes from 
the complex measure; Aq and Aq itself are always real-valued. 



2.3. Analytic continuation 

The sign problem of finite density QCD is manifested here as an imaginary term in the EQCD 
action, Eq. (|2,ip . This makes the standard Monte Carlo importance sampling impractical, 
except for very small chemical potentials and/or small volumes. One option to circumvent 
this problem is to use analytic continuation to complex values of Jl: the sign problem vanishes 
for purely imaginary p,. 

However, we emphasize that the direct analytic continuation in ft is clearly suboptimal and 
unnecessary in this case: of the terms appearing in EQCD action Eq. (j2T]) . only z7 3 TrL4g] 
is responsible for the sign problem. Thus, it is sufficient to analytically continue 73 (or z) to 
imaginary values and leave the other parameters to the values determined by the desired value 
of fi. By far the dominant effect of non-zero £1 is due to the /2-dependence of the parameter 
y in Eq. (|2.9p . we can take into account almost all of the effects of the chemical potential by 
just using the correct y(fi). The remaining small corrections are then taken into account by 
analytic continuation z — > iz0 

1 In [26] the susceptibility was evaluated by ignoring this correction; the improved statistics here make the 
small correction non-negligible. 
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Because the action (|2.ip is invariant under the simultaneous change z — ► — z and Aq — > — Ao, 
the partition function must be an even function of z (and fx). From this follows that the 
expectation values (-Aq) are even (odd) functions of z for even (odd) n. Therefore, we can 
Taylor expand the condensates Cj appearing in the expression for the susceptibility (I2.13h in 
powers of z as appropriate: 

d(z) = J2 Ci, n z n = i n Ci, n {-iz) n - (2.14) 

n n 

The analytic continuation now proceeds as follows: we perform simulations with imaginary 
value of z and determine the Taylor series coefficients n for each of the condensates up to 
the desired order. Using Eq. (|2. 14[) we obtain the the condensates Cj at real values of z, 
which can be inserted in Eq. (|2. 13|) in order to obtain the susceptibility. 

The dependence of the condensates on z is very mild, as expected, and it turns out to be 
sufficient to expand the condensates to very low order: 

Ci = ai + a 2 z 2 C 3 = a 4 

dC 

C 2 = a 3 C 4 = ^ = -2ia 2 z. (2.15) 

Note that we assume that C 2 and C3 are independent of z. This is indeed the case to the 
statistical accuracy we can reach. 

If we now denote with Ci{z\) the condensates measured from simulations with imaginary 
z = (0, zi), the susceptibility at real z = (zr,0) becomes 

X3,ij 



+ Trlkjij vl C 3 (zj) + + flj) y — C 4 {zi) . 

lO 47T z\ 



We note here that one simulation at some Z\ is sufficient to obtain the condensates and the 
susceptibility at all (small enough) zr. However, because both y and z depend on p,, only 
the value of zr, which corresponds to used in evaluating y is physical. Thus, for each value 
of the chemical potential we need to do a new simulation. We also choose to use z\ = in 
our simulations, eliminating the ratios zj/zr in Eq. (|2.16p . In what follows we shall use the 
notation z = zr = z\ to refer to both quantities. 

The phase diagram of EQCD has 3 distinct phases: a symmetric phase with (-Aq) = 
and 2 broken phases with non-zero (Aq), related by reflection (Aq) «-> — (Aq) |28j. In order 
to properly represent 4d QCD, the effective theory must remain in the symmetric phase. 
In the absence of the chemical potential the symmetric phase is at most metastable, when 
the parameters x and y are fixed to values which correspond to 4d QCD. This is normally 
not a problem, because the metastability is very strong and for all practical purposes the 
symmetric phase remains stable. 



5 



Applying imaginary chemical potential to the full action would decrease the value of the 
parameter y(fi), Eq. f|2.9[) . Hence, the metastability would be reduced and finally completely 
lost at some value of imaginary ^. However, for our method of analytic continuation this 
problem is completely avoided: because we calculate y(/j,) with real [i, the value of y increases 
as fi increases. Thus, the physical symmetric phase remains stable at all values of ji. 

2.4. Relation to 4d physics 

The relation between X3,ij and the physical 4d susceptibility is given by 

where Ap = pqqd — Peqcd is the perturbative 3d<-^4d matching coefficient for pressure. 
This is perturbatively computable order-by-order in coupling constant expansion, because all 
perturbatively problematic infrared singularities of high temperature QCD are fully contained 
in EQCD. The matching coefficient is currently known to order 0(g 5 ) |17]F1 

The simulation results in Sec. [3] indicate that the 0(g 6 ) and higher order contributions 
to the matching coefficient are very small; indeed, if we compare our results with the 4d 
simulation results, we obtain an excellent fit when we assume that these contributions vanish. 
Thus, the 0(g e ) and above contributions to the susceptibility are strongly dominated by the 
contributions coming from EQCD. 

Because EQCD is derived using perturbation theory, the final results depend on the per- 
turbative scale Aj^jg. We shall use here the value A^jg = 245 MeV, which has been obtained 
from lattice simulations with 2 light Wilson quarks [30]. For the critical temperature we use 
T c = 170 MeV, yielding the ratio T c /A m = 0.7§ The comparison between EQCD and 4d 
QCD simulation results is somewhat sensitive to the precise value of this ratio, but it can 
vary ±10% without significantly affecting the quality of the match. The value 0.7 turns out 
to be close to the optimal one for the matching. 

Due to the perturbative nature of the matching equations it turns out to be convenient to 
do the matching by subtracting the known 3d perturbative susceptibility and adding the 4d 
one: 

^2-^3^3 -X 3 ) + ^2~- (2.18) 

Here xlT ari( ^ X pcrt are 3d and 4d perturbative results. We also note that the quantities 
Xuu = Xdd and Xud = Xdu are related to those used in [S] by 



Xq = 2( Xuu +Xud) (2.19) 
1 
2 



XI = ^(Xuu-Xud) (2.20) 



5 4 . . 

XC = gXuu - -Xud- (2.21) 



For the pressure the matching coefficient has been calculated to 0(g ) in a much simpler theory in Ref. [29] 

MS 



3 We obtain the same value by using the results rgT c — 0.438 [3T] and foA^jg- = 0.62 [30] 
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2.5. On the lattice 



The theory in Eq. (12. ip is discretized in a standard way, as described in p£j. Due to the 
superrenormalizability of the 3d theory the couplings A3 and g% do not run, and m| has 
well-known linear and logarithmic divergences as the lattice spacing a —* 0. When these 
divergences are subtracted the continuum limit is straightforward. 

The evaluation of the quark number susceptibility requires the measurement of the con- 
densates in Eq. f)2. 13j) on the lattice. Due to the superrenormalizable nature of the theory, 
measurements can be rigorously converted to MS scheme in the lattice continuum limit; be- 
cause MS was used in in the perturbative matching to 4d QCD, this also allows us to compare 
to 4d results. 

The relations between the condensates on the lattice and in continuum can be written in 
the limit the lattice spacing a — > (or j3 = 6/(53 a) ~~ * 00 ) as C3 E2] 

Ci,MS = °iA ~ clP - c 2 (ln/3 + cY) + 0(1/ p), 
C 2 ,m = C2,a - [C2 (ln/3 + c' 2 )] + 0(1//?), 
C 3m = C 3 , a + 0(l/(3), 
C 4m = C^ a + 0(l/(3). 

Here labels MS and a indicate that the quantity is calculated in MS or lattice regularization, 
respectively. The numerical coefficients are 

ci ps 0.1684873399, 

c 2 = ps 0.1519817755, 

0.66796(1), (2.23) 



^2 



c 2 = — « 0.0316628698900405, 

167T 2 

c' 9 ps 0.08848010. 



3. Lattice simulations 

The lattice simulations were carried out using two massless quark flavors (iVf = 2). We used 
nine different values of temperature T, varying from T ps 1.9A^jg up to ~ 9 x lO^A^g. The 
temperature values are shown in Table [3l While the largest temperature is huge in physical 
units, in 3d parameters the variation is much milder; this is related to the fact that QCD 
approaches weakly coupled theory at high T extremely slowly. Thus, an extreme range of 
high temperatures is required in order to reliably assess the convergence to the perturbation 
theory. 

At each temperature we use 6 values for z = (fi u + ^d)/(3vr 2 T), /i u = /id, up to z = 0.15 or 
(J, u /T ~ 2.22. This amounts to 54 different (T,fi) pairs. For each physical point simulations 



7 



rp jrp 

1 /J-c 


„, ( r r\ 


x(l ) 


1.32 


0.357 


0.13 






n i n 
u. 1U 


11.5 


0.711 


0.06 


204 


1.18 


0.035 


3600 


2.02 


0.020 


2.4 x 10 7 


3.09 


0.013 


6.2 x 10 9 


3.99 


0.010 


1.9 x 10 13 


5.31 


0.0075 


6.1 x 10 16 


6.62 


0.006 



Table 1: The temperatures and corresponding i/q, x-values used in the simulations. For each 
temperature quark chemical potential has 6 values, parametrized by z = 2fi q / '{J>tt 2 T) = 
0,0.025,0.05,0.075,0.1,0.15, and for non-zero z y is modified according to Eq. (12. 9p . At 
each (T, /x)-pair the simulations are done using 6 different lattice spacings, parametrized by 
= 6/(sf o) = 32, 40, 54, 67, 80, 120. 

are done using 6 lattice spacings, parametrized by f3 = 6/ (<7§a) = 32 . . . 120. Thus, the lattice 
spacings vary by almost a factor of 4, enabling reliable extrapolation of the continuum limit. 
For the smallest lattice spacing (j3 = 120), the largest lattice size varies between 256 3 to 320 3 . 

In addition to the simulations at physical parameter values, we also did a several series of 
runs at fixed x, y and varying z. While these simulations do not correspond to any physical 
parameter set, they enable us to look at the z-dependence of the condensates separately. All 
in all, our dataset contains 693 individual runs. 

3.1. Continuum extrapolation 

It turns out that the accuracy requirement are so high that the continuum limit extrapolation 
of the condensates have to be taken with great care. Especially the continuum extrapolation of 
(-Ap) is critical, because it strongly dominates the susceptibility. While we know the divergent 
(as a — » 0) and constant contributions appearing in the continuum limit, Eq. (|2.22p . 0(a), 
O(alna) and higher order terms are not yet known. Thus, we use an ansatz 

(Al) a -C.T. = c 1 + C j + C jlo g (P) + ^, (3.1) 

where C.T. indicate the known counterterms in Eq. (|2.22p and Cj are fit parameters. The 
existence of the logarithmic term in the ansatz increases the errors of the extrapolation an 
order of magnitude compared to the case without the logarithmic term. However, c' 2 is 
expected to be a constant independent of y: by dimensional grounds the expansion of (Aq) 
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in powers of the lattice spacing can be written as 

{Al) a = ^ + D 2 gl + a [D 3 g 4 3 + D 4 m 2 + D 5 g 2 X + D 6 A 2 ] + 0(a 2 ). (3.2) 

The form of the 0(a) coefficients -D4, -D5 and Dq is known and they do not contain a term 
logarithmic in a |16|, I32|. whereas the coefficient D 3 , which is constant in y (or might 
include one. Thus, the possible a In a -contribution should indeed be independent of y. 

The existence of the logarithmic term can be seen in Fig. [H where we show the parameter c 2 
obtained from continuum fits using the ansatz (|3.ip . Note that here c' 2 is fitted independently 
for each physical parameter set, allowing arbitrary y (and z) dependence. As expected, the 
result is fairly well consistent with constant c 2 ~ 0.69; the remaining systematic discrepancies 
in the fit can be caused by contributions which are of higher order than 0(a 2 ), including 
terms of type a 2 In a. Thus, we shall fix c' 2 to this value in Eq. (|3.ip for all continuum limit 
extrapolations which follow. 

We note that the value of c' 2 has negligible effect on the results at small (physically relevant) 
temperatures; c 2 could be set to zero without affecting the continuum limit. It is significant 
only at very large T, where it potentially has a role when we compare simulations with the 
perturbation theory. We observe deviations from perturbative results even at very high T 
if c 2 < 0.4 (see also |26|). However, variations of order ~ 15% around 0.69 do not affect the 
final results. 

Nevertheless, it is clear that an analytic calculation of 0(a) effects in EQCD would be 
highly desirable. There is an ongoing calculation using stochastic perturbation theory |33| . 
which will hopefully confirm our results. 

The contributions of the other condensates are numerically much smaller and we were not 
able to see any sign of logarithmic a-dependence in those. It turns out that it is advantageous 
to make the continuum extrapolation using the full expression of the susceptibility f|2. 16|) . 
instead of extrapolating individual condensates. (Naturally, after the subtraction of the 
known counterterms in Eq. (|2.22p .) This extrapolation is shown in Fig. [2 

The lattice volumes are chosen large enough so that finite volume effects become negligible. 
We have tested this by doing simulations at selected parameter values using different volumes; 
at the smallest lattice spacing (J3 = 120) the volume varies from V = 144 3 up to V = 320 3 . 
No systematic finite volume effects inside two sigma errors. For more discussion of finite size 
effects on a related model see Ref. [2~T] . 

The finite chemical potential dependence is studied using the method described in Sec. 12.31 
The condensates C 2 = Vgl(((A 3 ) 2 ) - (i 3 ) 2 ) and C 3 = Vgl(((A 2 ) 2 ) - (A 2 ) 2 ) should be 
largely independent of z (for fixed x, y) for the equation (j2.15j) to be valid. This indeed 
turns out to be the case, within the statistical errors, and any remnant z-dependence is 
completely drowned out by the contributions from z-independent parts in Eq. (|2,15p . Indeed, 
the overall z-dependence of each of the condensates in Eq. (|2.15p turns out to be statistically 
almost invisible, with the exception of C4 = Vg 2 {(A^A^ — (Aq)(Aq)), which has a linear 
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Figure 1: Fitting of the logarithmic coefficient c' 2 in continuum extrapolation. The data is 
consistent with the assumption that c' 2 is a constant x 2 /d.o.f ~ 57/35 
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z-dependence. In practice the //-dependence of the susceptibility is almost completely due to 
the //-dependence of the parameter y and the /z 2 C3-term Eq. (|2.16|l R Nevertheless, here we 
do take into account the small z-dependence of C\ = {Aq) and C4, although it will affect the 
final results by only about 1 sigma. 

3.2. Diagonal susceptibility 

Now we are in position to compare the continuum limit results with the perturbation the- 
ory. First we shall look at the diagonal susceptibility X3,uu = X3,dd- The susceptibility has 
been calculated in perturbation theory up to order g e \nl/g [9]. In 3-dimensional units the 
perturbative result can be written as a power series in 1/^/yo, with the following result: 

XZI 8 + 9fi 2 3y 3 J 2 
gl V4 + 3fi 2 4vr 

(9 - 30x)/Z 2 - 4(3 + lOx) + 6(4 + 3/Z 2 ) ln(4 + 3/t 2 ) + 6(4 + 3fi 2 ) ln(y ) 3y 

2(4 + 3/i 2 ) (4^)2 

_ (8 + 3/Z 2 )(89 + 4tt 2 - 441n(2)) 9^ /2 
8(4 + 3fi 2 ) 3 / 2 (4tt) 3 

+ |576[-3438 + 40(2 + fi 2 )fi 2 ] + 119313vr 2 + 640(4 + 3fi 2 ) 2 ln(4 + 3fi 2 ) 

4- 640(4 4- 3P r mfa.)} 144(4 + 3 p 2)2(47r)4 4- 3^A. 2 + 0(y^\ (3.3) 

We have set here flu = fid = fi- As can be observed in Fig. [3] the overall agreement between 
the lattice result and the perturbation theory is very good, especially at large y (large tem- 
perature). The result contains an unknown //-independent order 0(?/q) -term denoted by /?M2 
in [18]. The same term appears also in the off-diagonal susceptibility, Eq. (|3.5p . and it turns 
out that it gives much tighter constraints for the value of /?M2 than the diagonal one. The fit 
to the off-diagonal susceptibility gives (3m2 = —0.1 ± 0.3. This value is small enough that its 
effect is negligible for the diagonal susceptibility, nonetheless we set here (3m2 = —0.1. 

In order to quantify the contributions not included in the perturbative result we calculate 
the difference x!j att — Xjj 61 * an d fit a function of form bi/y + 62/2/0 to the result. The fit 
results are shown in Fig. [J] and Table [2j We note that for small z the 1 /y -term is much 
smaller than the 1/yo-term, indicating that the 0(g 7 ) -contribution arising from EQCD is 
smaller in magnitude to the 0(g 8 ) term, at least for all physically relevant temperatures. At 
large z the statistical errors grow rapidly; this is due to the term oc z 2 V "^(((Aq) 2 ) — (A 2 ^ 2 ) 
in Eq. (I2TT6]) . 

Finally, we obtain the physical 4d result for the diagonal susceptibility from Eq. (|2.18p . 
As described in Sec. 12. 4^ the 3d<->4d mapping remains sensitive to the unknown 0(g & ) and 
higher order perturbative contributions to the matching coefficient. In Fig. [5] we show the 

4 This fact was used in the preliminary results published in ref. [27] 
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Figure 3: Left: the diagonal quark number susceptibility X3,uu/53 a t different values of 
chemical potential. The symbols indicate the lattice measurements, and the solid lines are 
the perturbative result. Right: The difference between the lattice and perturbation theory. 



z 


fit 


X 2 /dof 





0.0008(8)/^ -0.0137(8)/y 


12/6 


0.025 


0.0016(8)/^-0.0143(7)/y 


18/6 


0.05 


0.000(l)/ v ^-0.012(2)/y 


11/6 


0.075 


-0.003(2)/^ - 0.007(2)/y 


17/6 


0.1 


-0.008(4)/^ - 0.002(3)/y 


7.5/6 


0.15 


-0.012(6)/Vyo + 0.004(5)/y 


9.7/6 



Table 2: Fitting a function of form b\j ^fy^ + &2/2/0 to (xjfu U ~~ Xjt^uVflf- The smallest y 
(lowest temperature) points are left out of the fit. 
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-0.005 - 
-0.01 - 

-0.015 
-0.02 

-0.025 
-0.03 



z = 



- -0.025 



_l I I I L_ 



z = 0.05 



0.4 0.8 1.2 1. 

-1/2 

y 




0.4 0.8 1.2 1 

-1/2 

y 



0.4 0.8 1.2 Lf 

-1/2 

y 




0.4 0.8 1.2 1 

-1/2 

y 



0.4 0.8 1.2 1. 

-1/2 

y„ 



0.4 0.8 1.2 l.t 

-1/2 

y 



Figure 4: The diagonal susceptibility (xlfuu - X3 uidV^/flf as a function of 1/y/y with different 
values of the chemical potential. Solid line is a 1st order polynomial fit. The data at y ~ 
1.6 have been excluded from the fit. 



EQCD data at p, = with these unknown contributions set to zero. We observe that the 
result fits the 4d lattice simulations very well, clearly indicating that the magnitude of these 
contributions must be small, and in what follows we shall set them to zero. On the other 
hand, it should be noted that the difference between the purely perturbative result and EQCD 
simulation result is substantial at T < 10T C , as indicated by the two lines in Fig. [5j This is 
a clear indication that the contributions beyond the currently known perturbative ones have 
non-negligible effect at experimentally accessible temperatures. 

The //-dependence of the diagonal susceptibility is shown in Fig. [6j normalized to the 
Stefan-Boltzmann value 

xM=T 2 + \fi 2 . (3.4) 

7T 

We note that at temperatures above 100T C the deviation from the Stefan-Boltzmann law 
is independent of /i, but at lower T there is significant //-dependence. The /^-dependence 
matches very well the 4d lattice results by Allton et al. [5|, also shown in Fig. [6l 
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Figure 5: The diagonal susceptibility Xuu m 4d units at fi = 0. The data points indicate the 
full EQCD result. The continuous line is the result of the fit in Table[2j The dashed line shows 
the perturbative result alone, Eq. (|3.3p . using the same matching as in the EQCD result. 
The difference between these two curves indicates the magnitude of the non-perturbative 
contributions. The agreement with the 4d-lattice results of Gavai et al. [2] and Karsch et al. 
[I] is good. 




Figure 6: Left: The diagonal susceptibility at different /i, normalized to Stefan-Boltzmann 
law. Right: /^-dependence of the susceptibility compared with the 4d lattice results of Allton 
et al. [5]. 
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3.3. Off-diagonal susceptibility 



The perturbative result for the off-diagonal susceptibility in 3d units is 




4^/4 + 3fL 2 4^ 
27/x 2 y 



+ 



4 + 3/2 2 (4vr) 2 

27/x 2 (89 + 4tt 2 - 441og(2)) y\ 



+ 



8(4 + 3£ 2 ) 3 / 2 (4vr) 3 



+ 



3(2047168 - 1193137T 2 + 15360/Z 2 )/2 2 + 2560(4 + 3/i 2 ) 2 ln(4 + 3/2 2 )+ 



2560(4 + 3^) I°fa»)} 576(4 + 3p2)2(4 , )4 + + Ofe^ 2 ) 



(3.5) 



where (3m2 is the same unknown coefficient which appears in the diagonal susceptibility, 
Eq. (13. 3ft . In this case we can fit the value at z = 0, obtaining 



This value is small enough to have in practice negligible effect on the final results. Again the 
simulation data is very well described by the perturbation theory, Fig. [7J only at z = or at 
lowest temperatures can we observe deviations from perturbation theory. 

After matching to 4d, we obtain the result for off-diagonal susceptibility Xud ; shown in 
Fig. At T> 10T C the results match the perturbation theory very well, but at lower tem- 
peratures there are deviations: most significantly, at T = 1.32T C and \x = the simulation 
results clearly undershoot the perturbation theory. On the other hand, the 4d lattice results 
in [5] at \i = indicate small but non-zero value, which agrees well with perturbation the- 
ory. [9l [10] . This can be an indication that this point is already outside the validity range 
of EQCD; however, we also note that by increasing T c /Aj^g the EQCD results are brought 
closer to 4d lattice results [5J. The agreement with the perturbation theory and 4d lattice 
results is rather good already at T = 2.3T C . 

We also note that the physical value of the off-diagonal susceptibility is obtained in EQCD 
by a subtraction of two divergent as a — * terms; thus, as opposed to full 4d QCD simulation, 
there is no natural approximate symmetry which would force it to be small. Therefore, if 
EQCD starts to approach the limits of the validity, one can expect substantial deviations 
from physical results, as seen at T = 1.32T C . 

Nevertheless, the overall //-dependence of Xud is i n rough accordance with the 4d lattice 
results [5j already at T = 1.32T C , as shown on the right panel in Fig. [HI and at T = 2.3T C the 
agreement is already very good. The non-diagonal susceptibility is seen to behave quite well 
up to large values of fi/T ~ 2. 



Pm2 = -0.1 ±0.3. 



(3.6) 
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Figure 7: Left: the off-diagonal susceptibility X3,udA?3 m 3d units. Right: the difference 
between lattice and perturbative susceptibilities Odfud ~~ xludV^i) shown at 2 smallest fL. 
The statistical errors grow rapidly as fi increases. 
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Figure 8: Left: The off-diagonal susceptibility in 4d. At low temperatures we obtain signif- 
icantly different values from the perturbation theory (solid lines), but there is no deviation 
anymore at T = WT C . Right: /i-dependence of off-diagonal susceptibility compared with All- 
ton et al. [5] . The precision of results from [5] are probably not accurate enough to predict 
the behaviour at region fj,/T > 1. 
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4. Conclusions 

We have measured the quark number susceptibility of high temperature finite density QCD 
using lattice simulations of EQCD, an effective 3-dimensional theory of full 4d QCD. The 
very good match to the 4d lattice results with 2 light quark flavors at low temperatures and 
with the perturbation theory at high temperatures shows the wide range of applicability of 
the method. The diagonal susceptibility is seen to agree with 4d simulations by Allton et 
al. [5] even below 2T C , including the dependence on ji. On the other hand, we observe a 
substantial deviation from the known perturbative result up to temperatures ~ 20T C . The 
off-diagonal susceptibility is compatible with perturbation theory already at T> 10T C . The 
results also agree with the 4d simulations [5] except perhaps at lowest temperatures, T < 2T C . 

The results clearly indicate that EQCD is a viable method to obtain quantitatively signif- 
icant results of the hot QCD plasma down to T ~ 2T C . Equally significant is the observation 
that the currently known perturbative result alone deviates significantly from the correct 
result: while the perturbative result can be made to match the 4d lattice data by adjust- 
ing the still unknown (high perturbative order) matching coefficients, EQCD allows us to 
directly measure the differences between simulations and perturbative calculations without 
any scale or matching ambiguities. Thus, simulations of EQCD are exceptionally well suited 
for observing the convergence of the perturbation theory. It is worth noting that while the 
EQCD susceptibility also suffers from matching ambiguity, we obtain an excellent fit to 4d 
simulations by assuming these matching coefficients vanish, indicating that the contribution 
from these is necessarily very small. 

5. Acknowledgements 

We acknowledge useful discussion with K. Kajantie, M. Laine and A. Vuorinen. This work 
has been partly supported by the Magnus Ehrnrooth Foundation, a Marie Curie Fellowship 
for Early Stage Researchers Training, and the Academy of Finland, contract number 114371. 
KR also acknowledges partial support by the National Science Foundation under Grant No. 
PHY05-51164. Simulations have been carried out at the Finnish IT Center for Science (CSC). 

References 

[1] M. Asakawa, U. W. Heinz and B. Muller, Fluctuation probes of quark deconfinement, 
Phys. Rev. Lett. 85 (2000) 2072 [arXiv:hep-ph/0003169| ; 

S. Jeon and V. Koch, Charged particle ratio fluctuation as a signal for QGP, Phys. Rev. 
Lett. 85, 2076 (2000) |arXiv:hep-ph/0003168] . 

[2] R. V. Gavai, S. Gupta and P. Majumdar, Susceptibilities and screening masses in two 
flavor QCD, Phys. Rev. D 65 (2002) 054506 [hep-lat/01 10032] . 



17 



R. V. Gavai and S. Gupta, The continuum limit of quark number susceptibilities, Phys. 
Rev. D 65, 094515 (2002) [arXiv:hep-lat /0202006| . 

F. Karsch, S. Ejiri and K. Redlich, Hadronic fluctuations in the QGP, Nucl. Physifcon. 
A 774 (2006) 619 [hep-ph/05101~26] . 

C. R. Allton, S. Ejiri, S. J. Hands, O. Kaczmarek, F. Karsch, E. Laermann and 
C. Schmidt, The equation of state for two flavor QCD at non-zero chemical potential, 
Phys. Rev. D 68 (2003) 014507 [arXiv:hep-lat/0305007| . 

C. R. Allton et al, Thermodynamics of two flavor QCD to sixth order in quark chemical 
potential, Phys. Rev. D 71 (2005) 054508 |arXiv:hep-lat /0501030] . 

C. Bernard et al., QCD thermodynamics with 2+1 flavors at nonzero chemical potential, 
Phys. Rev. D 77 (2008) 014503 [ arXiv:0710. 13301 [hep-rat]]. 

C. Bernard et al. [MILC Collaboration], QCD thermodynamics with three flavors of 
improved staggered quarks, Phys. Rev. D 71 (2005) 034504 [arXiv:hep-lat /040 5029 1 . 

Y. Maezawa et al., Thermodynamics and heavy- quark free energies at finite temperature 
and density with two flavors of improved Wilson quarks, PoS LAT2007, 207 (2007) 
|arXiv:0710.0945l [hep-lat]]. 

M. Bluhm and B. Kampfer, Flavor Diagonal and Off-Diagonal Susceptibilities in a 
Quasiparticle Model of the Quark-Gluon Plasma, arXiv:0801.4147 [hep-ph]. 

A. Vuorinen, Quark number susceptibilities of hot QCD up to g**6 ln(g), Phys. Rev. D 
67 (2003) 074032 [hep-ph/02 12283] . 

J. P. Blaizot, E. Iancu and A. Rebhan, Quark number susceptibilities from HTL- 
resummed thermodynamics, Phys. Lett. B 523, 143 (2001) larXiv:hep-ph/0110369|. 

A. Ipp, K. Kajantie, A. Rebhan and A. Vuorinen, The pressure of deconfined QCD 
for all temperatures and quark chemical potentials, Phys. Rev. D 74, 045016 (2006) 
[arXiv:hep-ph/0604060"l . 

P. H. Ginsparg, First Order And Second Order Phase Transitions In Gauge Theories At 
Finite Temperature, Nucl. Phys. B 170 (1980) 388. 

T. Appelquist and R. D. Pisarski, High- Temperature Yang-Mills Theories And Three- 
Dimensional Quantum Chromodynamics, Phys. Rev. D 23 (1981) 2305. 

K. Kajantie, M. Laine, K. Rummukainen and M. E. Shaposhnikov, Generic rules for 
high temperature dimensional reduction and their application to the standard model, 



Nucl. Phys. B 458 (1996) 90 hep-ph/9508379]. 



E. Braaten and A. Nieto, Free Energy of QCD at High Temperature, Phys. Rev. D 53 
(1996) 3421 [hep-ph/9510408] . 

K. Kajantie, M. Laine, K. Rummukainen and M. E. Shaposhnikov, 3d SU(N) + 
adjoint Higgs theory and finite-temperature QCD, Nucl. Phys. B 503 (1997) 357 
[arXiv:hep-ph/9704416| . 



18 



K. Kajantie, M. Laine, K. Rummukainen and Y. Schroder, The pressure of hot QCD 
up to g 6 \n(l/g), Phys. Rev. D 67 (2003) 105008 [hep-ph/021132l] . 

A. Vuorinen, The pressure of QCD at finite temperatures and chemical potentials, Phys. 
Rev. D 68 (2003) 054017 [hep-ph/03051~83 |. A. Vuorinen, The pressure of QCD at finite 



temperature and quark number density, arXiv:hep-ph/0402242, 



A. Gynther and M. Vepsalainen, Pressure of the standard model at high temperatures, 
JHEP 0601, 060 (2006) [hep-ph/0510"37 5|. 

K. Kajantie, M. Laine, K. Rummukainen and Y. Schroder, How to resum long- 
distance contributions to the QCD pressure? Phys. Rev. Lett. 86 (2001) 10 
[arXiv:hep-ph/0007109| . 

A. Hietanen, K. Kajantie, M. Laine, K. Rummukainen and Y. Schroder, Plaquette 
expectation value and gluon condensate in three dimensions, JHEP 0501 (2005) 013 
[hep-lat/0412008] . 

A. Hietanen and A. Kurkela, Plaquette expectation value and lattice free energy of three- 



dimensional SU(N(c) gauge theory, hep-lat/0609015 



M. Laine and Y. Schroder, Two-loop QCD gauge coupling at high temperatures, JHEP 
0503 (2005) 067 | arXiv:hep-ph/0503061] . 

K. Kajantie, M. Laine, J. Peisa, A. Rajantie, K. Rummukainen and M. E. Shaposh- 
nikov, Non-perturbative Debye mass in finite T QCD, Phys. Rev. Lett. 79 (1997) 3130 
[arXiv:hep-ph/9708207| . 

A. Hart, M. Laine and O. Philipsen, Static correlation lengths in QCD at high tempera- 
tures and finite densities, Nucl. Phys. B 586 (2000) 443 [hep-ph/0004060| . 

A. Hart, M. Laine and O. Philipsen, Testing imaginary vs. real chemical potential in 
finite-temperature QCD, Phys. Lett. B 505, 141 (2001) [arXiv:hep-lat /0010008| . 

A. Hietanen and K. Rummukainen, Quark number susceptibility at high temperature, 
PoS LAT2006 (2006) 137 [arXiv:hep-lat/0610111| . 



A. Hietanen and K. Rummukainen, Quark number susceptibility of high temperature and 
finite density QCD, PoS LAT2007, 192 (2007) [ arXiv:0710.5 058 [hep-lat]]. 

K. Kajantie, M. Laine, A. Rajantie, K. Rummukainen and M. Tsypin, The phase 
diagram of three-dimensional SU(3) + adjoint Higgs theory, JHEP 9811 (1998) 011 
[arXiv:hep-lat/9811004] . 

A. Gynther, M. Laine, Y. Schroder, C. Torrero and A. Vuorinen, Four-loop pressure of 
massless O(N) scalar field theory, JHEP 0704 (2007) 094 [arXiv:hep-ph/070330 7 1 . 

M. Delia Morte, R. Frezzotti, J. Heitger, J. Rolf, R. Sommer and U. Wolff [ALPHA 
Collaboration], Computation of the strong coupling in QCD with two dynamical flavours, 
Nucl. Phys. B 713 (2005) 378 [arXiv:hep-lat/0411025| . 



19 



[31] V. G. Bornyakov, S. M. Morozov, Y. Nakamura, M. I. Polikarpov, G. Schierholz and 
T. Suzuki [DIK Collaboration], Finite temperature Lattice QCD with two flavors of im- 
proved Wilson fermions, PoS LAT2007 (2007) 171 [ arXiv: 071 1.14271 [hep-lat]]. 

[32] M. Laine and A. Rajantie, Lattice- continuum relations for 3d SU(N)+Higgs theories, 
Nucl. Phys. B 513, 471 (1998) | arXiv:hep-lat/9705003] . 

[33] C. Torrero, M. Laine, Y. Schroder, F. Di Renzo and V. Miccio, Towards J^-loop NSPT 
result for a 3-dimensional condensate- contribution to hot QCD pressure, arXiv:0711. 1176l 
[hep-lat] . 



20 



